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Abstract - Methods for determining the percolation threshold usually study the behavior of 
network ensembles and are often restricted to a particular type of probabilistic node/link removal 
strategy. We propose a network-specific method to determine the connectivity of nodes below 
the percolation threshold and offer an estimate to the percolation threshold in networks with 
bidirectional links. Our analysis does not require the assumption that a network belongs to a 
specific ensemble and can at the same time easily handle arbitrary removal strategies (previously 
an open problem for undirected networks). In validating our analysis, we find that it predicts the 
effects of many known complex structures (e.g., degree correlations) and may be used to study 
both probabilistic and deterministic attacks. 



Introduction. — The study of percolation in com- 
plex networks has broad applications including epidemic 
spreading [l] , propagation of excitation in neural networks 
[2], and robustness of networks to random failure [3] or 
strategic attack [ij[5j . A central problem is estimating the 
percolation threshold, the critical fraction of nodes or links 
of an initially connected network that must be removed to 
disintegrate it into small disconnected fragments. Knowl- 
edge of how a network fragments can improve strategics 
for designing attack [4|5] and immunization techniques [6] 
or increasing network robustness [Tj[8J . 

Several studies have proposed techniques to estimate the 
percolation threshold of a network for various situations 



4U9 - 13 . These studies typically use ensemble approaches, 



where one studies the typical behavior of a set of networks 
satisfying some set of properties. Common ensembles in- 
clude: (i) networks with fixed degree sequence generated 
from, for example, the configuration model [14] ; (ii) net- 
works with an expected degree sequence generated with, 
for example, the Chung-Lu model 15 ; and (iii) Marko- 



vian networks with correlations between nearest neigh- 
bors, where correlations may be captured by the prob- 
ability P(d'|d) that a node with degree d is connected 
to a node of degree d'. Here d = (d m ,d out ) denotes the 
the number of incoming (d m ) and outgoing (d out ) links at 
a given node. While significant progress has been made 



in the study of such ensembles [4{[9}|T2], we note two im- 
portant limitations of the ensemble approach, (i) Given a 
single network, it is not clear what ensemble should be se- 
lected to capture the properties of the network. Typically, 
networks found in applications contain various structural 
properties (e.g., correlations 12 , clustering 16 , and com- 



munity structure 17 ) that are not always accounted for 



in the ensembles. A related problem is that it has been re- 



cently observed 18 that, given an ensemble of networks, 



some network properties can vary significantly within the 
ensemble. Thus, it is not clear that ensemble approaches 
give the best description of a single network found in prac- 
tice, (ii) Some ensemble theories are impractical when ap- 
plied to individual real networks. For example, theories 
to estimate the percolation threshold in Markovian en- 
sembles |11| require the estimation of a potentially dense 
and very large, (CO X (Ca.C'.L) matrix P(d|d'), 
which is difficult to obtain from a single network. 

In this Letter, we develop a network-specific theory 
for percolation in networks with bidirectional links. In 
contrast to previous approaches to network-specific esti- 
mates of the percolation threshold which were restricted 
to purely directed networks 



13 , our work has a broader 



range of applicability as it can be applied to undirected 
networks and directed networks with bidirectional links. 
Our method is based on an analysis of the network's 
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adjacency matrix, which is often known or can be es- 
timated in important applications (e.g., the power grid 
19 and air transportation networks [20])- Besides relax- 
ing the ensemble assumptions of previous research (e.g., 
that the network is strictly Markovian) , one significant 
advantage of this approach is that it can easily account 
for arbitrary strategies of node/link removal. Network- 
specific approaches are therefore well suited for developing 
network-specific attack/defense strategies, immunization 
techniques, etc. In addition to estimating the percolation 
threshold, we predict the expected number of nodes ac- 
cessible to each node after the network disintegrates. This 
has various applications such as predicting the outbreak 
size of an epidemic [l]. We finally show that our method 
may be used to study the fragmentation of a network sub- 
ject to either probabilistic or deterministic attack. 

Analysis. — We formalize weighted percolation (i.e., 
in which nodes and/or links are retained with arbitrary 
probabilities) as follows: for a network with N nodes 
described by a possibly asymmetric adjacency matrix A 
(A nm = 1 if a link exists from node n to node m and 
A nm — otherwise), node n is retained with probability 
q n , and the directed link from node n to node m is retained 
with probability p nm . Letting q = N^ 1 J2 n In denote the 
average node retention probability, qN nodes are expected 
to remain after a realization of this process (referred to 
as a percolation trial). Unweighted node percolation cor- 
responds to q n = q < 1 and p nm = 1, while unweighted 
link percolation corresponds to q n = 1 and p nm = p < 1. 
For our analysis, it is useful to introduce a matrix A with 
entries defined by qA nm — q m PnmA nm , which represents 
the probability that a link exists from node n to node m, 
given that node n is retained. Because our analysis de- 
pends only on the matrix A, it is applicable to link, node, 
and mixed (i.e., simultaneous link and node) percolation. 
However, for the remainder of this Letter we consider only 
node percolation, p nm = 1. 

For a given node-targeting strategy, defined as a set of 
retention probabilities {q n }, we are interested in the size 
s of the largest strongly- connected component (LSCC), the 
largest subset of nodes so that any node in the subset is 
reachable from any other node in the subset. The perco- 
lation threshold q* is defined as the value of q such that 
s -C N for q < q* (the subcritical regime) and s ~ N for 
q > q* (the supercritical regime). For a particular removal 
strategy, the subcritical regime may be analyzed by noting 
that after a percolation trial, only a fraction of the network 
is reachable from a given node n following directed links. 



Following 10 we define, for a given percolation trial, the 



out-component of node n as the set of nodes that may be 
reached from node n via the remaining network (including 
node n) and define s°"* as the size of the out-component 
of node n averaged over many percolation trials. To moti- 
vate subsequent analysis, consider first the case when the 
network is a directed tree. In that case s° u * satisfies the re- 



counts the nodes reachable from node n by counting the 
nodes reachable from its neighbors, and adds 1 to account 
for node n itself. The same expression approximately ap- 
plies to directed networks that are locally tree-like 21 

mH where A is 



13 



and leads to the results in rcf. 
the principal, or Perron-Frobenius, eigenvalue of A (i.e., 
Au = Au). When links are allowed to be bidirectional, 
however, the expression above overestimates the size of 
s°"', since the terms s^' on the right hand side might 
include nodes that are reachable by following links back 
into node n (see fig. fll). 




Fig. 1: When computing sj" , to compensate for the over- 
counting of nodes due to the bidirectional link 1 O 3, /3i3 
reduces the contribution of s 3 ut on s° ut . We approximately 
have fa ~ 3/4 and sf*' = 1 + qA 12 s° 2 ut + qA 13 fi 13 s 3 mt . 

To correct for this over-counting of nodes, we heuris- 
tically modify the contribution of s^f on the right hand 
side by a factor /3 nm (to be determined), 



nmMnrn^rn 



(i) 



To determine a self-consistent expression for (3 nm , we 
note that from eq. the relative contribution of the 



nmHnm' J rn 



■hi 



out-component of node m on s° is qA 
Therefore, to reduce the contribution of s™ 1 on s° ut to 
account for the branch returning to node n (if present), 

j™'. Inserting he 
and solving for /3. 



we let p nm = 1 - qA mn /3 mn s° ut / s™* . Inserting here the 
corresponding expression for j3. 
obtain 



we 



Pnm (1 q A nm A rnn ) 



1 - 



qA, 



After substitution of eq. ^ into eq. ([!]) , we find 
s out = [I-D(q)]- 1 y, 



(2) 



(3) 



where s out = [sf 1 
N, y is a vector with entries 

.2 \ ^ 



s^ 1 '] , / is an identity matrix of size 



Vn 



1 



A nk A kn (l-q 2 A nk A kn )- 1 
and D(q) is a matrix with entries 



k 



D nm (q) = qA nm y n (l 



a 2 A A 



(4) 



(5) 



Given a removal strategy, eq. ^ can be solved to obtain 
the expected out-component size for each node. To obtain 
an estimate for the percolation threshold, note that eq. Q 
requires the invertibility of the matrix / — D(q). This 
matrix is invertible when Xntq) < 1, where XdM is the 



lation s°"* = 1 + X)m 94mC f i where the right hand side principal eigenvalue of D(q). As \D( q ) — > 1 the out- 
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component sizes diverge as s° Mt ~ [1 — Xjj( q - ) ]~ 1 w n , where 
w is the principal eigenvector of D(q). A similar argument 
can be made for the divergence of the in-component sizes. 
Since the LSCC above the percolation threshold can be 
thought of as the set of vertices with infinite in- and out- 
components [11], we predict the percolation threshold as 



q * D = '■ Xd ^ = 

qe[0,l] 



(6) 



We note that if there are no bidirectional links, 
A nm A mn — and D(q) — qA, and the results of ref. [13] 
are recovered. While one may solve eq. ^ numerically, it 
is both practical and insightful to approximate eqs. ^ 
[6| for large s° u * and small q. Letting s° Mt 3> 1 and 
finm ~ 1 in eq. yields the approximate eigenvalue prob- 
lem s° ut w qj2 m inmC*- lt follows that s out cx u. Upon 
substitution we find q ~ A^ 1 under these conditions, yield- 
: 1 - A- 1 ^ 



ing to first order /3„ 



lUn/um. Defining 



n m ■ A nrn I 1 



a,, 



\u r , 



(7) 



with principal eigenvalue equation Cz = Acz and using 
y « 1 = [1, 1, . . . , 1] T , we obtain the predictions 

-li 



(/ - q cy 



q c 



Ac 1 - 



(8) 
(9) 



In addition to offering simplified predictions for s° ut and 
q* , for unweighted percolation (i.e., A = A and u = u) 
these estimates allow us to bound Ac using the prin- 
cipal eigenvalue A of the network adjacency matrix A 
(e.g., Au = Au). Direct application of the Bauer-Fike 
Theorem 24 for the limiting case of an undirected net- 
work yields | Ac - A| < WX^UAU' 1 ^ = 1, where U = 
diag[ui, . . . , ujv]- Finally, considering l T Cu and using 
z ~ u yields Ac ~ A — 1. One implication of these results 
is that q* — > for large A, which is consistent with the lack 
of an unweighted percolation threshold for well-connected 
networks such as scale-free networks |9j. We note that 
eqs. Q and ([£]) are in best agreement with eqs. ^ and 
^ near q = q* and when the network is strictly undi- 
rected or strictly directed. 

Examples. — In what follows, we will motivate the 
need for our theory, recover previous results, explore sev- 
eral applications, and illustrate the robustness of our anal- 
ysis to complex structures in networks. We will consider 
both computer-generated and real-world networks. 

We first highlight the need for a network-specific 
method for undirected networks and show that unlike the 
ensemble approach, a network-specific method captures 
variability in q* across an ensemble. We consider per- 
colation in an uncorrelated, random network formed by 
retaining the giant component from an Erdos-Renyi pro- 
cess [22j (henceforth referred to as an ER network) with 
N = 10 4 nodes and 37V links. In fig. |] we show the frac- 
tion of retained nodes in the LSCC, s/N, as a function 
of q (solid grey curve). Our predictions for q* given by 
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Fig. 2: (Colour online) The relative size of the giant compo- 
nent, s/N (solid grey line), is shown for unweighted percolation 
on an ER network |22| (see text). Our predictions to q* given 



by eq. |6| and eq. |§| are given by vertical solid black and 
dashed red lines and are in good agreement with the undirected 
ensemble result |10| (green squares') . T he network-specific pre- 
diction for directed networks, A -1 [li] (blue dot-dashed line), 
is shown to be inaccurate. The inset shows experimental values 
for q* vs. eq. Q for networks within an ensemble. Unlike the 
ensemble approach, variation in q* is naturally accounted for 
by our network-specific approach. 



eq. ffify (black line) and eq. d9| (red dashed line) and the 
undirected ensemble theory [10] (green squares) work well, 
whereas the network-specific theory for directed networks, 



A 
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(blue dot-dashed line), does not, as expected. The 
inset shows experimental values for q* versus our predic- 
tion using eq. ([6| for an ensemble of uncorrelated networks 
obtained by rewiring the ER network while retaining a 
fixed degree sequence, [di, g?2, • • • , djv]. The rewiring pro- 
cess is similar to that in refs. [5]|T2], except an additional 
step is taken to ensure the resulting network contains all 
nodes in its LSCC. Note that whereas our network-specific 
method naturally accounts for variability in q* across the 
ensemble, the ensemble approach (squares) cannot, pre- 
dicting q* rj 0.33 for all members of the ensemble. 

Experimental values for q* in this inset and remaining 
figures were found using extrapolation over 6 overlapping 
intervals of width 0.03 that span [0.03,0.13] to find the 
intercept with the horizontal axis. This yielded a mean 
and standard deviation (shown when significant). While 
the actual percolation threshold q* is often well defined 
for model networks in the N — > 00 limit (e.g., by using 
traditional finite-size scaling arguments), one can only es- 
timate it for real-world networks which have a fixed and 
finite size. We also note that while q* can often be ex- 
tracted by examining the size of the second-largest cluster, 
this approach was observed to significantly overestimate q* 
for the relatively small networks considered here for which 
N < 10 5 . 

We next illustrate our analysis in networks with bidi- 
rectional links and recover previous results for unweighted 
percolation in uncorrelated networks. Letting / denote 
the fraction of directed links, we begin with an undirected 
(i.e., f = 0) ER network with N = 10 4 nodes and 5N 
links and iteratively replace randomly chosen undirected 
links with directed links in a random orientation until the 
network is strictly directed (i.e., f = 1). One can ob- 
serve in fig. [3] that our predictions by eqs. (|6| and ^ 
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Word-Association Network 




(d(d-l)} 



Fig. 3: (Colour online) Equations |6| and eq. (JoJ agree well 
with ref. 'Ill and experimental values for predicting the un- 
weighted percolation threshold q* of an uncorrelated random 
network. / denotes the fraction of links that are directed, which 
increases from (undirected) to 1 (directed). 



agree well with the ensemble result 11 and experimental 
values. Note that the predictions of eq. Q agree with 
those of eq. ([6| very well for f — or 1 and that max- 
imal disagreement occurs near / ~ 0.7. We note that 
one disadvantage of extrapolating on a fixed range (e.g., 
s/N £ [0.03,0.13]) while varying q* is that a drift occurs 
due to the combination of finite-size affects and the fact 
that the slope of an s(q) curve just above q* varies mono- 
tonically with q* . This is one cause for the disagreement 
between theory and experiment in fig. [3] for large /. 

We now recover the ensemble predictions for strictly 
undirected and strictly directed networks shown by hori- 
zontal lines in fig. [3j For / = 0, our result Ac ~ A — 1 and 
the undirected mean-field (MF) result A « (d 2 ) / (d) lead to 
q* « (d)/(d(d-l)) (9). For / = 1 we have A nm A mn = 
and C = A, which recovers q* sa A -1 
MF result A -1 « (d m )/(d m d out 
of ref. 
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ran 

Again, the 
13 recovers the result 



13 



Our next examples explore targeted attacks, where we 
let q n oc d l n for I € {-1,0, 1}. Here d n = «™ + d°^)/2 
denotes the degree of node n averaged over incoming and 
outgoing links. Preferentially removing nodes with large 
degrees (i.e., I < 0) can model probabilistic attacks [2] and 
biased infrastructure failure 3 , whereas preferentially re- 
moving nodes with small degrees (I > 0) may represent 
a non-obtrusive degradation. We consider real-world net- 
works for which previous methods are either impractical 
or not justified for the reasons mentioned in the introduc- 
tion. In fig. [4ji we show the relative size of the LSCC, 
4> = s /(qN), f° r weighted percolation on a directed Word- 
Association (WA) Network 25 . Vertical dashed lines and 
x's represent the predictions of eq. §6§ and eq. The 
inset shows this graph with the horizontal axis respec- 
tively normalized by the value of q* found from eq. ([6]) 
for each targeting strategy, where a LSCC of size s ~ N 
appears at q/q* ~ 1 for all curves. In fig. |4j3, we show 
the prediction by eq. ^ (solid lines) and observed (sym- 
bols) values of q* as a function of I, normalized by the 
value of q* at I = 0, for the WA network 25 and an 
undirected network of Facebook (FB) friendships at Cal- 
tech 



20 



Being a well-connected network with power-law 
degree distribution, the FB network has a very small un- 
weighted percolation threshold [q*(Q) ~ 0.01, not shown]. 
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Fig. 4: (Colour online) (a) Average values for cj> = s/(qN) over 
32 trials are shown for weighted percolation with q n oc d l n on 
the scale-free WA network 25 . Vertical dashed lines and x's 



represent the predictions of eq. |6| and eq. Q. The inset shows 
the same quantities with the g-axis normalized by eq. |6| for 
each I value, (b) The dependency of q* on / is shown for the 
WA (circles) and FB networks 26 (x's). 

However, preferentially removing nodes with large degree 
can yield nontrivial thresholds with nearly a ten-fold in- 
crease [q*(-l) ~ 0.1] . 

Our next example illustrates eqs. ^ and ((8| in a peer- 
to-peer (P2P) network of file downloads [23] with N = 
6301 nodes for weighted site-percolation with q n oc g?~ 5 . 
In this example s° u * can be interpreted as the expected 
number of nodes infected by a computer virus released by 
user n and preferentially targeting nodes with large degree 
[l] . In fig. [5] we show a sample of the values of s° u * pre- 
dicted by eq. Q (x's) [which agrees with eq. ([8|] versus 
their experimental values, s° u * (exp). The prediction given 
by eq. ^ (x's) is very accurate for the average removal 
probability q = 0.2 shown in fig. [5^,. For the larger value 
q = 0.4, our prediction deviates somewhat from the ob- 
served values as shown in fig.[5]D. This is expected because 
s° u * is predicted to diverge at q* ~ 0.476, but experimen- 
tal values are bounded by the finite network size V, so 
the predicted value must become larger than the observed 
value as q — > q* . 
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Fig. 5: ( Colo ur online) (a) For weighted percolation on a P2P 
network [23] with retention probability q„ oc d~ ' 5 and average 
retention rate q = 0.2, the values of s£ u * predicted by eq. (|3| 
(x's) agree well with experimental values, s™* (exp), which 
were averaged over 2 16 percolation trials, (b) Similar results 
for q = 0.4, where some deviations appear (see text). 
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Fig. 6: Malicious attacks were simulated on two significant 
undirected networks: (a) a power grid for the U.S. 
flight network for US Airlines 



19 and (b) a 



20 , where the node with largest 



dynamical importance [5] was iteratively removed (along with 
nodes disconnected from the LSCC after its removal). The rel- 
ative size of the LSCC is shown in the top panels (grey lines) 
as a function of the fraction of removed nodes, k/N. The rel- 
evant eigenvalues for the networks' adjacency matrices, i.e., 
A (dashed grey lines) and Ac (solid grey lines) are shown in 
the bottom panel as a function of the fraction k/N of removed 
nodes. Black lines indicate the same variables under random 
node removal. The disappearance of the LSCC corresponds to 
Ac ~ 1 (horizontal dotted lines). 



In our next example we apply our results to network 
fragmentation under deterministic attack. We consider 
two important undirected networks, a power grid [19| and 
an airline transit network |20| , and point out that knowl- 
edge of how they fragment is essential for taking measures 
of protection. Simulated attacks were implemented by it- 
eratively removing the node n corresponding to the largest 
dynamical importance, DI n oc u n v n [5 j, where u, v are the 
right and left principal eigenvectors of A. For these undi- 
rected networks, DI„ = v? n . After each node removal, com- 
ponents fragmented from the LSCC were also removed. In 
the upper panels of fig. [6j we show the relative size of the 
LSCC after the removal of k nodes (grey lines), whereas in 
the lower panels we show both A (dashed grey lines) and 
A c (solid grey lines), where C is given by eq. 0. The 
horizontal dotted line denotes Ac = 1. For comparison, 
black lines indicate the the same variables for random node 
removal. The disappearance of the LSCC corresponds to 
Ac ~ 1. For these undirected networks this corresponds 
to A — 1 ~ 1, whereas for directed networks we recover the 
result of ref. [13] , A ss 1. Here A is the principal eigenvalue 
of the network adjacency matrix A. 

In the remaining examples we demonstrate the robust- 
ness of our results to complex structures in networks. Be- 
ginning with degree correlations, we first address Marko- 
vian correlations between the degrees of adjacent nodes, 
which can be characterized by the assortativity coefficient 
r e (- 



-1,1 12 



where r > (r < 0) indicates that nodes 
tend to connect to other nodes with similar (different) de- 
grees. We again allow the retention probability of a node 
to depend on its degree, q n cx d l n . We note that ref. [4] pro- 




Fig. 7: (Colour online) (a) Values of q* predicted from eq. |6| 
(solid lines) and observed experimentally (symbols) for perco- 
lation on an ER network having TV = 10 4 nodes and 5iV undi- 
rected links rewired with correlations (measured by assortativ- 
ity coefficient r). The three curves correspond to unweighted 
percolation (I — 0, middle), and weighted percolation (I = — 1, 
top, and 1 = 1, bottom). Squares show the undirected Marko- 
vian prediction |12| . (b) Eq. (J©J) (solid line) predicts observed 
values of q* (stars) for a non-Markovian network (see text), 
whereas the Markovian prediction (squares) cannot. 



vides an ensemble approach for a similar removal strategy, 
but their analysis is restricted only to scale-free networks 
lacking correlations. In fig. [7^, we show the effect of de- 
gree correlations on q* for an ER network with N = 10 4 
nodes and 5 AT undirected links that is rewired to have 
correlations. Experimental values are shown for I = — 1 
(x's), I = (stars), and I — 1 (crosses). Solid lines indi- 
cate the prediction of eq. which was found to coincide 
with that of eq. ([9|, and the squares indicate the undi- 
rected Markovian prediction of ref. [12] (applicable only 
for 1 = 0). Degree correlations were varied while keeping 
the degree distribution constant following the algorithm 
in refs. 5, 12 . Note that while assortativity promotes ro- 
bustness for unweighted percolation (/ = 0) by reducing 
q* , we find that its effect is reduced (amplified) for I < 
(I > 0). For example, the percolation threshold is largely 
unaffected by degree corralations for q n cx d^ 1 for this 
network (see x's in fig.JT^i). 

Turning to degree-correlations of the non-Markovian 
type in undirected networks, a case which no previous the- 
ories can handle, we consider the set of networks used to 
produce fig. [7^, but subject them to the following rewiring 
process: each link m«nis replaced by two new links and 
a new node j, m«j« n, resulting in correlations across 
paths of length two. In fig. [TJa, whose horizontal axis is 
carried over from fig. [7^,, we show that the network-specific 
prediction eq. ^ (solid line) agrees with the observed val- 
ues of q* (stars) for unweighted percolation on these non- 
Markovian networks. For comparison, direct application 
of the Markovian ensemble method [12] (squares) does not 
give good results. 

We conclude by showing that our results remain accu- 



rate for relatively large clustering coefficient, c 16 . In this 



experiment links were iteratively added to an ER network 
by finding paths of length two and completing the trian- 
gles. In fig. [8]we show our predictions for and the observed 
value of q* as a function of c. 
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Fig. 8: (Colour online) Predictions to q* and experimental val- 
ues for percolation on an ER network with N = 10 4 nodes 
and 5iV links as a function of the clustering coefficient, c (see 
text). Note that the lack of degree correlations yields good 
agreement between our results, eq. (6) and eq. (9), and the 
ensemble result of ref. 
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Conclusions. — We have presented a network-specific 
approach to weighted percolation in undirected networks 
and directed networks with bidirectional links (previously 
open problems) . As opposed to most previous theory deal- 
ing with network ensembles, our method predicts unique 
percolation characteristics for each unique network. While 
ensemble and network-specific methods offer complemen- 
tary strategies which may lead to different insights, in this 
Letter we highlighted several benefits of a network-specific 
approach, (i) Ensemble methods cannot account for vari- 
ability across networks within the ensemble (see fig. [2] and 
ref. [18]). (ii) Application of any ensemble approach to a 
real-world network requires a priori assumptions about 
that particular network (e.g., that it either lacks complex- 
ity not accounted for in the ensemble or that its effects are 
small), (iii) Our approach naturally accounts for degree 
correlations of the non-Markovian type, (iv) Arbitrary 
targeting for node and/or link removal can be easily han- 
dled with our approach, (v) A network-specific analysis 
allows one to study deterministic attacks. Our results help 
explain why strategically decreasing and increasing A of- 
fer fundamental strategies to respectively attack [5] and 
protect [7] networks. 

We have provided several examples showing that our 
results are robust to various attack strategies, degree dis- 
tributions, degree correlations, and moderate clustering. 
However, when strong community structure is added we 
expect our analysis to breakdown as occurs with the en- 
semble approaches [27] . Finding a network-specific theory 
accounting for communities remains open for future re- 
search. Another direction of future work includes applying 
our techniques to more complicated percolation problems 
such as k-core percolation and Achlioptas processes. How- 
ever, it is likely that efforts following the techniques pro- 
vided here will also be restricted to the subcritical regime. 
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